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Evidence of the Dynamics of Relativistic Jet Launching in Quasars 

Brian Punsly^ 

ABSTRACT 

Hubble Space Telescope (HST) spectra of the extreme ultraviolet (EUV), the 
optically thick emission from the innermost accretion flow onto the central su- 
permassive black hole, indicate that RLQs tend to be EUV weak compared to 
the radio quiet quasars (RQQs); yet the remainder of the optically thick ther¬ 
mal continuum is indistinguishable. The dehcit of EUV emission in RLQs has a 
straightforward interpretation as a missing or a suppressed innermost region of 
local energy dissipation in the accretion flow. This article is an examination of 
the evidence for a distribution of magnetic flux tubes in the innermost accretion 
flow that results in magnetically arrested accretion (MAA) and creates the EUV 
dehcit. These same hux tubes and possibly the interior magnetic hux that they 
encircle are the source of the jet power as well. In the MAA scenario, islands of 
large scale magnetic vertical hux perforate the innermost accretion how of RLQs. 

The hrst prediction of the theory that is supported by the HST data is that the 
strength of the (large scale poloidal magnetic helds) jets in the MAA region is 
regulated by the ram pressure of the accretion how in the quasar environment. 

The second prediction that is supported by the HST data is that the rotating 
magnetic islands remove energy from the accretion how as a Poynting hux dom¬ 
inated jet in proportion to the square of the fraction of the EUV emitting gas 
that is displaced by these islands. 

Subject headings: Black hole physics — magnetohydrodynamics (MHD) — galax¬ 
ies: jets—galaxies: active — accretion, accretion disks 


1. Introduction 

The mechanism that drives powerful beams of radio emitting plasma, moving near the 
speed of light, in radio loud quasars (RLQs) has been the subject of much speculation in the 


^1415 Granvia Altamira, Palos Verdes Estates CA, USA 90274 and ICRANet, Piazza della Repubblica 
10 Pescara 65100, Italy, brian.punslyl@verizon.net 



2 


literature, with little or nothing in the way of observations of the region that launches the 
jets (Lovelace 1976; Blandford and Znajek 1977 Blandford and Payne 1982 [Punsly 2008). 
This circumstance has now changed based on Hubble Space Telescope (HST) observations 
of the extreme ultraviolet (EUV) in quasars ( Punsly||2014 ). Quasars are generally associated 
with the optically thick thermal emission from gas that accretes onto a supermassive black 


hole 

Lynden-Bell and Rees 

(1971); 

Shakura and Sunyaev 

(1973) 

Novikov and Thorne 

(1976); 

Malkan 

(1983) 

Szuszkiewicz et al. 

(1996 

). Curiously, ~ 10% of these quasars have conspic- 


uous beams, or jets, of relativistic plasma on scales that can exceed a million light years and 
powers 10^ - 10"^ the integrated light of the largest galaxies, RLQs. The optical/ultraviolet 
spectra of RLQs and radio quiet quasars (RQQs) tend to be very similar except for subtle 


differences in certain emission line strengths and widths (Steidel and Sargent 1991 Boro- 


son and Green 1992 Brotherton et ah 1994). These emission line regions are far from the 


central engine, ~ 10^ - 10^ larger than the central black hole radius. Consequently, this 
research path has provided very little understanding of the jet launching mechanism. The 
EUV continuum, A < 1100 A, is created orders of magnitude closer to the central engine 


and RLQs display a significant EUV continuum deficit relative to RQQs (Telfer 2002). 


The quasar luminosity is widely believed to arise from the viscous dissipation of turbu¬ 


lence driven by the differential rotational shearing of accreting gas (Shakura and Sunyaev 


1973). In numerical and theoretical models, the highest frequency optically thick thermal 


emission arises from the innermost region of the accretion flow and its frequency is shortward 


of the peak of the spectral energy distribution (SED) (Zhu et ah 2012). Consider this in 
the context of the RLQ and RQQ quasar composite spectra from HST data in Figure 1 


(Telfer 2002). The continuum of the composite spectra are indistinguishable except for the 
EUV emission shortward of the peak of the SED at ~ A = 1100A (where the spectra are 
normalized to 0). Thus, the difference in the EUV emission between RLQs and RQQs likely 
arises from suppressed emission in the innermost region of the accretion flow in RLQs of 


what is otherwise a similar accretion flow to that found in RQQs (Punsly 2014) 


In this article, the explicit predictions of a MAA (magnetically arrested accretion) de¬ 
scription of the EUV deficit in RLQs are explored (see Figure 2 for a schematic illustration). 
It is posited that islands of large scale magnetic vertical flux perforate the accretion flow of 


RLQs within a few black hole radii of the central black hole Punsly (2014); Igumenshchev 


et al. (2003); Igumenshchev (2008). Three pieces of information are synthesized to constrain 


the dynamics in the inner accretion flow, the long term time average jet power, Q, the peak 
luminosity of the SED (a good surrogate for the bolometric luminosity, Lboi, of fh® accretion 
flow) and the EUV dehcit. MAA indicates that the rotating magnetic islands remove energy 
from the accretion flow as a Poynting flux dominated jet in proportion to the square of the 
fraction of the EUV emitting gas that is displaced by these islands and in direct proportion 
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to the accretion rate. These predicted relationships are well fit by the data extracted from 
the HST spectra, thus providing strong support for the MAA interpretation of RLQs and the 
paradigm that magnetic flux in the inner accretion flow is the switch that launches quasar 
jets (|Meier||1999|). 


In the following an extended magnetically arrested scenario is explored (EMAA) in 
parallel with the analysis of the MAA scenario. In the EMAA model, the magnetic flux that 
threads the EUV region is the outermost extent of an interior magnetic flux distribution 
that threads the equatorial plane of the ergosphere and the event horizon of the central 
supermassive black hole. The one necessary assumption that is required to connect this 
interior region to observation is that the vertical magnetic flux in the interior region is 
proportional to the magnetic flux in the EUV region, oc ^euv- Iii this circumstance, the 
interior magnetic flux distribution is a major source of jet power. The EUV deficit is created 
by the flux that perforates the innermost accretion flow. Based on the assumption above, 
this interaction also provides an indirect tracer (a probe with a linear response function) of 
With this assumption, the observational data discussed here also supports the EMAA 
model. 


The paper is structured as follows. The second section is a review of MAA and EUV 
suppression. In this section, the fundamental relationships between jet power and EUV 
suppression are found. The goal of this article is to test these predicted relationships versus 
observations. In order to define the experiment, a method for estimating the jet power 
from the radio data and a sample of quasars with suitable HST spectra and radio images 
need to be established. These are the subjects of Sections 3 and 4. Section 5 describes the 
experimental tests performed and the conclusions. The last section is a discussion of the 
results. 


2. Jet Production and Suppressed EUV 


The dynamics of viscous shear in black hole accretion flows with and without MAA 
are reviewed in order to assess the affects of MAA in an annular ring of the accretion flow 
surrounding the central black hole. Assume that magnetic islands fill a fraction, fy, of the 
volume of the ring, V, and penetrate a fraction, /, of the top and bottom surface areas 
of the annular volume, SA. The volume of magnetic islands is Vmi and its complement 
in V is V^j, fv = J (IVmi/V. The surface area elements of the top and bottom faces 
are, cISAmi and dSA'^j, respectively. The turbulent dissipation that heats the plasma in 
accretion flows is produced as a consequence of the magneto-rotational instability (MRI) in 


the 3-D numerical simulations De Villiers et ah (2003); Penna et al (2010). The accretion 
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Composite Quasar SED 



Log Frequency (Hz) 


Fig. 1.— The EUV Deficit. The HST quasar composite continuum accretion disk spectra. 
The blue and red plots are the ultraviolet SED peak and the EUV composite spectra for 
RQQs and RLQs, respectively (Telfer et al 2002). Note the steeper SED (larger oeuv) for 
RLQs. 
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Fig. 2.— Magnetically Arrested Accretion. The logarithmic false color contonr plot of the 
vertical poloidal magnetic pressnre (i.e., Bp = Bz) from an MAA simnlation is nsed to 
schematically illnstrate the concepts posited in this stndy. One can see Fignre 2 of Pnnsly 
et al (2009) for the color scale, bnt this is not necessary for this discnssion. The false 
color accentnates the two-flnid natnre of the accretion flow. The red to yellow regions are 
magnetic islands, the origin of strong Poynting jets and weak EUV emission. The dark bine 
(the densest plasma) and green are strong EUV emitting regions. 
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flow in MAA simulations is perforated by large scale magnetic flux tubes, magnetic islands 
(see Figure 2). Thus, MAA creates a two component fluid. Firstly, there are regions in 
which energy is removed from the fluid by turbulent dissipation. Secondly, in the magnetic 


islands, energy is removed from the how vertically as Poynting hux ( 

Igumenshchev et ah 

2003 

Igumenshchev||2008 

McKinney et al.||2012; 

Tchekhovskoy et ah 

2011 

). The large scale 


magnetic flux suppresses the MRI induced dissipation in these regions. Simulated MAA flows 


are subsonic and therefore do not produce signihcant gas heating from shocks (McKinney et 


al. 

2012; 

Punsly 

2014 


Thus, viscous dissipation is the primary source of heat creation at 
the boundary of the magnetic islands and in surrounding accreting gas. The total volume 
available for MRI induced viscous heating is reduced by the magnetic islands. MAA provides 
an alternative to local turbulent dissipation for removing energy from the plasma allowing 
the gas to accrete closer to the black hole. The magnetic islands radiate Poynting flux, S^, 


along the magnetic held lines at the expense of the energy of the plasma (Igumenshchev 


2008). The local physics that produces the turbulent viscosity, rji, in V^j is unchanged from 


standard accretion and therefore so is the viscous stress in the huid surrounding the islands, 

( |2014| ), in the magnetically arrested case, the 


= rjtr^dTt/dr). It was shown in Punsly 


radiative luminosity is ~ 1 — / of what it wou 


d be for standard accretion with the same mass 


accretion rate. If / is the fraction of the inner accretion how surface area, SA, penetrated by 
magnetic islands in MAA then the EUV luminosity, L(EUV), obeys the approximate scaling 


Punsly (2014) 


L(EUV) oc (1 - f)SA . 


( 1 ) 


For any MHD Poynting hux dominated jet, regardless of the source, the total integrated 
electromagnetic poloidal energy hux is 


P 


( 2 ) 


where <I> is the total magnetic hux enclosed within the jet, dA^ is the cross-sectional area 
element and fc is a geometrical factor that equals 1 for a uniform highly collimated jet 


(Punsly 2008). Thus, not only do the magnetic islands of large scale poloidal hux in the 


inner accretion how suppress radiation from this region, but they provide a source of Poynting 
hux (power for the jet) as they orbit around the black hole with an angular velocity, flp. If 
M is the black hole mass in geometrized units, Vlp ~ c/M for rotation based on black hole 
spin or Keplerian orbits. Thus, Equation (2) implies for magnetic islands 


Q 


S^dA^ oc Bif FA oc ^A oc f^M{Vr/B) , 


( 3 ) 


where is the vertical magnetic held at the disk surface. The string of proportionality 
statements in Equation (3) require elaboration. The hrst proportionality arises from the 
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fact that magnetic flux, $ is conserved in the perfect magnetohydrodynamic approxima¬ 
tion. So the value in the jet is the same as that in the MAA region. In the MAA region 
$ = J BzdSAMi- This integral is approximated as $ ~ fBzSA. At this point it is worth 
reiterating the intent of exploring the EMAA model as well. Equations (2) and (3) can be 
evaluated using $ ^int for the jet power, and oc ^euv ~ fBzSA. Thus, the following 
analysis can be used to compare the pure MAA scenario or the EMAA scenario with ob¬ 
servation. The only qualihcation is that the EMAA analysis has the additional assumption 
that ^int scales linearly on average with ^euv- 


The second proportionality derives from pressure balance at the interface of the mag¬ 
netic islands and the enveloping accretion flow. In order to evaluate the pressure balance, 
it is important to realize that the magnetic islands are not frozen into the enveloping ac¬ 
cretion flow. For example, in Figure 2 the magnetic islands are in the process of spiraling 
outward, slowly. The dynamics are not time stationary. Outward motion of the stronger 
magnetic islands begins in the innermost regions of the accretion flow where the magnetic 
islands tend to merge (Igumenshchev p008 ). Furthermore, the gas density in the magnetic 
islands decreases due to outflow (i.e., a jet) in the inner accretion flow. The density is < 1% 
of the enveloping accretion flow ( [Igumenshchev 2008). Hence, they are weakly affected by 
gravity. By contrast, the dense enveloping accretion flow is strongly attracted to the central 
black hole. This is a classic Kruskal- Schwarzschild instability, the magnetized version of 
a Rayleigh-Taylor instability (Stix 1992). Thus, the magnetic islands become buoyant and 


drift outward relative to the radial ram pressure of the inflow as discussed in Punsly et 


al. (2009[) and the associated video of a simulated flow. The magnetic held in the islands 


decreases (expansion and ablation) as they drift outward in response to the decreasing dy¬ 
namic pressure in an attempt to achieve a new pressure balance with the enveloping medium 
(Igumenshchev 2008). This dynamic is unsuccessful and the islands slowly drift outward. 
In the MAA model considered here, the magnitude of the outward migration velocity, v^j^j 
is much smaller than the magnitude of the bulk radial velocity of the inward enveloping ac¬ 
cretion how, Vr- This is generally true in the simulations of Igumenshchev (2008), except at 
the smallest radii where the stronger islands are sometimes particularly unstable and burst 
outward until a more stable equilibrium is achieved. By contrast, the vertical magnetic hux 
that permeates the inner accretion how in the simulation of the Kerr geometry discussed in 
detail in Punsly et al. (2009); Hawley and Krolik (2006| seems to persist at the smallest radii 


for longer than it does in the simulations of Igumenshchev (2008), remaining within r < 2M 
(in geometrized units), the ergosphere, for ~ 2 local Keplerain orbital periods Punsly (2007). 
Designate quantities in the magnetic island by the label MI and in the enveloping accretion 
how by AF. There are four pressure components, Pg, Pb, Pr, and Pram, corresponding to 
gas, magnetic, radiation and ram pressure, respectively. Pr is a slowly varying quantity 















































around horizontal closed loops in a local neighborhood of the optically thick accretion flow 
that encircles each magnetic island. The radiation pressure is continuous through the low 
density magnetic islands. Thus, from the second moment of the radiative transfer equation, 
one expects Pr{MI) ^ Pr{AF) at the interface (Punsly 1996). In the inner accretion flow, 
Tram Pg ln luminous Quasars. Unlike Pr, the ram pressure is not continuous across the 
boundary, P^sj^{MP) P^i,^{AF), because of the large density differential and the motion of 
the magnetic islands relative to the enveloping accretion flow ( Punsly|1996 ). There is also no 
buoyancy force in the azimuthal direction. The force on the magnetic islands imposed by the 
inflow of matter is the strongest force resisting the buoyant outflow of the magnetic islands. 
Thus, despite all the uncertainty in the precise physics of magnetic island time evolution (see 
the hnal section), the radial ram pressure should determine the internal magnetic pressure 
of the islands to hrst approximation. Therefore, the pressure balance at the interface of the 
magnetic island and the enveloping accretion flow is approximately P^axaiAP) k. Pb{MI) 
(Igumenshchev 2008). Consequently, the substitution oc Pram in Equation (2) was made. 
In summary, if the magnetic islands are not extremely short-lived transient features in the 
inner accretion flow, an approximate balance of the ram pressure and the magnetic pressure 
of the poloidal magnetic held in the islands must exist. 


The third proportionality on the right hand side of Equation (2) arises from mass con¬ 
servation near the black hole, where the “pseudo-half angle”, 9, is the ratio of cross sectional 
area to SA (evaluated near the outer boundary of the MAA region). Since, as discussed 
above, Vr — v\jj ~ to the accuracy of scaling laws in Equation (3), Pram cc v^. Consider 
this in the context of the bolometric luminosity of the accretion disk, Lboi = 'r]{a)MF, where 
MF is the mass-energy accretion rate and the accretion efficiency is r]{a), a function of black 
hole spin, a. Equation (3) can be transformed into an approximate relationship that is more 
conducive to comparison with observation 


Q/F,i^C{l/v{a)){vr/0)f , 


(4) 


where C is a constant. For simplicity, it is assumed that the vast majority of black holes in 
quasars are rapidly rotating (Bardeen 1970 Elvis, Risaliti and Zamorani||2002 ). Secondly, it 
is also assumed that variations in 6 and Vr provide random scatter to Equation (4) and are 
not the primary physical drivers of Q. It is implicit in the approximate form of Equation 
(4) that the exact distributions of magnetic islands and the associated functional variation 
of UiT’ in the magnetically arrested region near the black hole are higher order corrections 
and only contribute cosmic scatter to the simplihed relationship. These approximations are 
collectively referred to as the homogeneous approximation for the MAA. This makes the 
following analysis tractable and not dependent of specihc models of the various parameters. 
With these assumptions. Equation (4) and the scaling with / of L(EUV) in Equation (1) 














implies an approximate relationship, 


Q/Lboi ^ A(1 - L(EUV)/L(EUV)rqq)2 , (5) 

where A is a constant and L(EUV)rqq is the hdncial EUV Inminosity if there were no 
magnetic islands. 


3. Long Term Time Averaged Jet Power. 


The two most viable options for estimating the jet power, Q, of qnasars are either based 
on the low freqnency (151 MHz) flnx from the radio lobes on 100 kpc scales or models of the 
broadband Doppler boosted synchrotron and inverse Compton radiation spectra associated 
with the relativistic parsec scale jet. Each method has its advantage. The 151 MHz method 
is generally considered more reliable since it does not involve large nncertainties dne to 


Doppler beaming (Willott et al 1999). A disadvantage is that it involves long term time 


averages, Q, that do not necessarily reflect the cnrrent state of qnasar activity. In this stndy 
Q = Q. A method that allows one to convert 151 MHz flnx densities, F 151 (measnred in Jy), 
into estimates of long term time averaged jet power, Q, (measnred in ergs/s) is captnred by 
the formnla derived in Willott et al (1999); Pnnsly| ( |2005 ): 

Q ^ [(f/15)3/2]i 1 X 10^® ergs/s , 

Z = 3.31 - (3.65) X 

^4 n ono \^3 ■ r\ Ar\ . r\ a a a rxr^r-l —0.125 


( 6 ) 


- 0.203X3 + 0.749X2 + 0.444X + 0.205] ’ 


(7) 


where X = 1+z, F 151 is the total optically thin flnx density from the lobes (i.e., contribntions 
from Doppler boosted jets or radio cores are removed). This sophisticated calcnlation of the 
jet kinetic Inminosity incorporates deviations from the overly simplihed minimnm energy 
estimates into a mnltiplicative factor f that represents the small departnres from minimnm 
energy, geometric effects, hlling factors, protonic contribntions and low freqnency cntoff 


Willott et al 

( !!)■)') 

). The qnantity, f, was farther determined to most likely be in the range 

of 10 to 20 

Blnndell and Rawlings 

(2000 

). In this paper we adopt the following cosmological 


parameters: Hq=70 km/s/Mpc, Da = 0.7 and = 0.3. Dehne the radio spectral index, 
a, as Fi, oc z/“". The formnla is most accnrate for large classical donble radio sonrces, thns 
we do not consider sonrces with a linear size of less than 20 kpc which are constrained by 
the ambient pressnre of the host galaxy. Alternatively, one can also nse the independently 


derived isotropic estimator in which the lobe energy is primarily inertial in form Pnnsly 


(2005) 
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Due to Doppler boosting on kpc scales, core dominated sources with a very bright one sided 
jet (such as 3C 279 and most blazars) must be treated with care (Punsly 1995). The best 
estimate is to take the lobe flux density on the counter-jet side and multiply this value by 2 
(bilateral symmetry assumption) and use this estimate for the flux density in Equations (6) 
-( 8 ). 


4. Sample Selection 


Determination of quasar EUV continua requires space based observations of modest 
redshift quasars since ground based observations of high redshift objects in which the EUV 


is redshifted into the optical are heavily attenuated by the Lya forrest (Zheng et ah 1997 


Telfer||2002 )■ In order to get a meaningful estimate of oeuv, a range of at least 700 A to 1100 
A in the quasar rest frame is needed to extract the continuum from the numerous broad 


emission lines (Teller 2002). Therefore, a redshift of > 0.63 is required. Troughs from 
Lyman limit systems (LLS) were removed by assuming a single cloud with a opacity. 
This was considered acceptable if the power law above the LLS could be continued smoothly 
through the corrected region (see the spectra in the Appendix). If there were many strong 
absorption systems or a LLS that compromised a broad emission line, this simple procedure 
was deemed inadequate for continuum extraction with the available data and the spectrum 
was eliminated from the sample. A small correction for the Lyman valley was also made 


(Zheng et ah 1997). Additionally, if there was evidence of a blazar synchrotron component 


contribution to the continuum such as a dominant flat spectrum radio core accompanied by 
high optical polarization or optical/UV variability, or low equivalent width of the emission 
lines, the underlying accretion disk continuum was considered too uncertain for the sample. 


As discussed in the last section, the most reliable methods of estimating the long term 
time averaged jet power Q are based on the the optically thin emission from relaxed radio 
lobes. Thus, all sources in the sample needed proof of extended emission on scales larger than 
the host galaxy so that the lobes can relax (> 20 kpc). Verification required archival high 
resolution interferometry images made between 0.408 GHz and 5 GHz. The HST and radio 
selection criteria resulted in a total of 20 sources for the sample. Note that two new southern 
hemisphere sources have been added to the sample of spectra from Punsly (2014). The Notes 
on Individual Sources in the Appendix provides the details that allow these sources to pass 
the criteria required to be in the sample. The optically thin emission was estimated based on 
151 MHz - 178 MHz flux densities (if available) and the lobe fluxes from the radio images. 
The largest spread in the estimates of Q, based on optically thin extended emission, are 


















Table 1: The Jet Power of the HST RLQ Sample 


Source 

Alias 

z 

Lobe Flux Density 
(mJy) 

Observed Frequency 
(GHz) 

Estimated Q 
(10^® ergs/s) 

0024+224 


1.12 

11 

1.4 

1.38 + 0.67 

0232-042 

PKS 0232-04 

1.45 

6100 

0.178 

32.45 ± 8.81 

0637 -7516 

PKS 0637 -752 

0.65 

209 

4.8 

2.83 + 1.24 

0743-67 

PKS 0743-67 

1.51 

1450 

2.5 

65.87+ 14.4 

0959+6827 


0.77 

27 

1.4 

0.27 + 0.16 

1022+194 

4C +19.34 

0.83 

2400 

0.178 

4.01 + 1.67 

1040+123 

3C 245 

1.03 

1179 

1.4 

15.18 + 4.95 

1137+660 

3C 263 

0.65 

18380 

0.151 

13.17 + 4.42 

1229-021 

4C -02.55 

1.05 

9000 

0.160 

20.52 + 6.25 

1241+176 

PG 1241+176 

1.27 

112 

1.4 

2.90 + 1.27 

1244+324 

4C +32.41 

0.95 

3370 

0.151 

6.53 + 2.51 

1252+119 

PKS 1252+119 

0.87 

15 

1.4 

0.31 + 0.18 

1317+5203 

4C+57.21 

1.06 

3800 

0.178 

10.38 + 3.66 

1340+289 

FBQS J1343+2844 

0.91 

830 

0.151 

1.67 + 0.79 

1340+606 

3C 288.1 

0.96 

9900 

0.151 

17.83 + 5.61 

1354+19 

PKS 1354+19 

0.72 

545 

1.4 

3.57 + 1.32 

1415+172 

PKS 1415+172 

0.82 

1020 

0.408 

3.31 + 1.42 

1857+566 

4C +56.28 

1.60 

8650 

0.160 

48.99 + 11.82 

2149+212 

4C +21.59 

1.54 

5300 

0.160 

30.20 + 8.35 

2340-036 

PKS 2340-036 

0.90 

61 

1.4 

0.78 + 0.41 


Table 2: The EUV Properties of the HST RLQ Sample 

Source 

Alias 

z 

\Lx(\ = llOOA) 
(10^^ ergs/s) 

^bol 

(10“^^ ergs/s) 

o^EUV 

Spectrograph/ 

Grating 

0024+224 


1.12 

64,05 

243.50 

1.77 ± 0.09 

FOS/G160L 

0232-042 

PKS 0232-04 

1.45 

94.47 

358.98 

1.75 ± 0.09 

FOS/G160L, G270H 

0637 -7516 

PKS 0637 -752 

0.65 

20.80 

79.02 

1.75 ± 0.09 

FOS/G160L 

0743-67 

PKS 0743-67 

1.51 

110.78 

420.95 

2.20 ± 0.30 

FOS/G190H, G270H 

0959+6827 


0.77 

22.28 

84.65 

1.40 ± 0.10 

FOS/G160L 

1022+194 

4C +19.34 

0.83 

7.02 

26.67 

2.55 ± 0.15 

FOS/G160L 

1040+123 

3C 245 

1.03 

11.90 

45.21 

2.20 ± 0.20 

FOS/G160L 

1137+660 

3C 263 

0.65 

27.39 

104.10 

2.00 ± 0.10 

COS/G130M, G160M; FOS/G190H 

1229-021 

4C -02.55 

1.05 

25.09 

95.36 

2.65 ± 0.15 

FOS/G160L, G190H, G270H 

1241 + 176 

PG 1241 + 176 

1.27 

86.78 

329.76 

1.79 ± 0.09 

STIS/G230L 

1244+324 

4C +32.41 

0.95 

13.77 

52.31 

2.41 ± 0.25 

FOS/G160L 

1252+119 

PKS 1252+119 

0.87 

16.12 

61.27 

1.45 ± 0.10 

FOS/G160L, G190H 

1317+5203 

4C+57.21 

1.06 

37.43 

142.25 

2.18 ± 0.20 

FOS/G160L 

1340+289 

FBQS J1343+2844 

0.91 

14.76 

56.08 

190 ± 0.10 

FOS/G160L 

1340+606 

3C 288.1 

0.96 

11.39 

43.27 

1.78 ± 0.09 

STIS/G140L, G230L 

1354+19 

PKS 1354+19 

0.72 

27.46 

104.34 

1.95 ± 0.10 

FOS/G160L 

1415+172 

PKS 1415+172 

0.82 

7.28 

27.66 

1.93 ± 0.10 

FOS/G160L 

1857+566 

4C +56.28 

1.60 

22.60 

85.87 

2.87 ± 0.14 

STIS/G230L 

2149+212 

4C +21.59 

1.54 

18.22 

69.24 

2.27 ± 0.20 

STIS/G230L 

2340-036 

PKS 2340-036 

0.90 

38.78 

147.35 

1.88 ± 0.09 

FOS/G160L 
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bounded on the high side by Equation (6) for the parameter, f = 20, and on the low side by 
Equation (8). These two extremes are used to generate the uncertainty in Q in Table 1. 


5. Experiment Design and Results 


Equations (4) and (5) are two approximate basic and testable predictions of the MAA 
explanation for the EUV dehcit in RLQs. In order to design an experiment to test the 
validity of the relations, one needs a sample of RLQs with adequate EUV spectra and radio 
observations to perform the test as described in the last section. Equations (4) and (5) 
are tested by looking at the HST spectra that cover the span from the SED peak at 1100 
A to 700 A. The first experiment tests Equation (4) by assessing the implied correlation 
between Q/Lboi and oeuv, where the spectral luminosity of the continuum is approximated 
as oc z/“"Euv_ Since the data used here covers the peak of the SED at A ~ llOOA, an 
accurate expression, Lboi ~ 3.8 ALa(A = llOOA), can be used to estimate the accretion disk 


luminosity proper (less reprocessed IR emission in distant molecular clouds) (Davis and Laor 
2011 Punsly 2014). The relevant data is captured in Tables 1 and 2. The data scatter is 


plotted in Figure 3. If the MAA scenario is correct not only would there be a correlation 
between cteuv with Q/Lbob but the correlation of cteuv with Q/Lboi should be stronger 
than ctEuv with Q since dividing by Lboi is equivalent to dividing out the scatter induced 
by the term in Equation (3). The Spearman rank correlation test indicates that the 
probability that the correlation of oeuv and Q/L^oi in Figure 3 occurs by random chance 
is 0.001. Note this is a more significant correlation compared to a probability of the scatter 
occurring by random chance of 0.006 for oeuv and Q and 0.150 for oeuv and redshift (z). 
Thus, the expected correlation exists. Physically, the stronger correlation of Q/Tboi with 
ttEuv compared to cteuv and Q is evidence that the poloidal magnetic pressure in the islands 
is being set by the ram pressure of the surrounding accretion flow. The correlation of Q/Tboi 
with oeuv is dehnitely improved, but not dramatically improved from that with Q and 
ctEuv- This is likely a consequence of the circumstance that the correlation with Q is already 
very strong. Allowing for cosmic scatter generating effects (i.e., a simplified uniform MAA 
assumption, modest geometric variations from object to object and epoch to epoch within 
the same object, plus non-contemporaneous measurement of Q with L{EUV)), it is probably 
not realistic to expect that the probability that the correlation of Q/Tboi with oeuv occurs 
by random chance is 0 to more than 2 significant digits. In order to see if this is a statistically 
signihcant result consider the partial correlation of Q/Lboi with oeuv when Q is held hxed. 
The partial correlation coefficient is 0.492 which correspondence to a statistical signihcance 
of 0.984. Conversely, the partial correlation of Q with oeuv when Q/L]^oi is held fixed is 
signihcant at the 0.581 level. Similarly, repeating this analysis with the Kendall tan rank 
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EUV Spectral Index Versus Jet Efficiency 

3.25 r 



-4.0 -3.5 -3.0 -2.5 -2.0 -1.5 -1.0 -0.5 0.0 


log (Q/Lboi) 

Fig. 3.— The EUV Deficit versus the Jet Efficiency. A scatter plot of Q/L^o\ and oeuv- 
The blue rectangle represents the RQQs for comparison. The RQQs composite of Figure 1 
has ctEuv = 1-57 ± 0.17 (Teller et al 2002). Another composite of predominantly RQQs has 
ttEuv = 1-41 ± 0.16 (Stevans et al 2014). So a broad range of 1.50 ± 0.25 is chosen to span 
both composites. 
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Normalized Jet Power Versus Normalized EUV Deficit 


0 I 



LOG(f) 


Fig. 4.— Jet Power and the EUV Deficit. The plot of log[Q/Lboi], normalized jet power, 
versns normalized EUV dehcit, log(/). The best £t line is plotted in black for the case that 
the lower cntoff on / is 0.025. The slope is 2.005 in agreement with the exponent of 2 in 
Eqnation (11). The determination of / for small valnes (flat EUV spectrnm sonrces) is very 
sensitive to the choice of the hdncial radio qniet EUV level, hence the large error bars for 
small /. 
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test, yields: the partial correlation of Q/L\,o\ with cteuv when Q is held hxed is statistically 
signihcant at the 0.957 level and the partial correlation of Q with oeuv when Q/Lboi is 
held fixed is statistically significant at the 0.753 level. There is statistical evidence that the 
correlation of Q/Tboi with oeuv is the fnndamental correlation and the correlation of Q with 
oeuv is a spurious correlation. This supports the premise of equation (3) that it is the jet 
power divided by accretion rate (ram pressure) that is the physical variable that is related 
to the EUV deficit. Finally, the much weaker correlation with redshift indicates that more 
profound physics is occurring, not just cosmic evolution or selection effects. 


The second experiment is to test Equation (5) which is much more restrictive than 
Equation (4) and much more sensitive to any initial assumptions and deviations from the 
homogeneous approximation. In particular, this experiment will test both the Q/Lboi de¬ 
pendence on p and the connection between / and the EUV deficit in Equation (1). First, 
define L(EUV) in a normalized form since from Equation (1) it scales with SA. The short¬ 
est wavelength that is uniformly sampled is A = 700A, so L(EUV) oc Ly{\ = 700A) is the 
best available proxy for L(EUV). In order to remove the dependence on SA this value is 
normalized by the peak of the SED 


L(EUV) 


L^(A = 700A) 

L^(A = llOOA) ■ 


(9) 


For RQQs, the average value of cteuv = 1-57 ± 0.17 from [Telfer (2002). As noted in Figure 
3, another composite of predominantly intermediate redshift RQQs has cteuv = 1-41 ± 0.16 
(Stevans et ah 2014). Thus, there is no unique estimate of the hducial or baseline EUV 
level for the absence of magnetic islands. The variation in the RQQ EUV translates into 
an uncertainty in the estimate for / in Equation (5). Denote the maximum and minimum 
estimates of the EUV in the absence of magnetic islands determined from the HST composite 
spectral analysis as 


L(EUV)rqq 


max 


L(EUV)rqq 


min 


= 700A) 
= llOOA) 
L^{X = 700A) 
L^{X = llOOA) 


a=1.25 


0=1.74 


0.57 , 
0.46 . 


( 10 ) 


Combining Equations (5), (9) and (10) yields the crude, but simple prediction of the homo¬ 
geneous MAA model applied to these intermediate redshift quasars observed with HST, 


Q/Lboi ^Af^A 


1 

■ L;.(A = 700A) ■ 

'L^{X = iiooA)' 

1 ^ 


.L^(A = llOOA). 

. L^(A = 700A) . 

RQQ 


( 11 ) 


The larger the estimate for the RQQ EUV baseline level in Equation (10), the larger the 
value of / (the EUV deficit) that is computed in Equation (11). Equation (11) is utilized as 
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follows in Figure 4 for estimating /: the expression for / in Equation (11) is computed two 
ways, one with the maximum baseline EUV luminosity from Equation (10) and again with 
the minimum baseline luminosity in Equation (10). The results of these two values in the 
expression for / in Equation (11) are averaged, this is the / value that is plotted in Figure 
4. The uncertainty in this average / is calculated as the maximum of the the expression for 
/ in Equation (11) minus the average. Implementing Equation (10) in Equation (11) is not 
completely straightforward since according to Table 2, two of the sources have oeuv <1-74 
and two others quasars have cteuv that are close to this value. The volume of the EUV 
emitting region that is displaced by magnetic islands cannot be negative, / cannot be less 
than 0 in Equation (11). Thus, consider the example of a lower bound or cutoff of at least 
2.5% of the accretion surface area in the EUV region (/ > 0.025) is displaced by magnetic 
islands for the RLQs with flatter EUV spectrum when the upper limit value of cteuv = 1-74 
is used in Equation (11) as the condition for the absence of magnetic islands. The result is 
shown graphically in Figure 4 for each source in Tables 1 and 2. This cutoff is arbitrary, so 
the effect of its variation is investigated below. 

The second experiment in Figure 4 plots the logarithm of Q/Tboi versus logarithm of / 
from Equation (11) from the data in Tables 1 and 2. The black line is the linear £t to these 
20 sources and the slope is 2.01. The theoretical prediction of MAA is 2 in Equation (11). 
The incredible agreement is clearly coincidental since there is significant scatter and many 
crude approximations in the derivation of Equation (11). As an illustration of this claim, 
consider a lower bound of / > 0.05. In this case, the slope is 2.22. Alternatively, if the lower 
bound of / > 0.01 is used, the slope is 1.88. The fit to the best fit line in Figure 4 is fairly 
good with a squared multiple regression correlation coefficient of 0.57. 


6. Implications and Conclusion 

The primary results of this study that are presented in Tables 1 and 2 (as well as Figures 
3 and 4) are empirical. It was shown that HST EUV spectra and radio interferometry data 
are consistent with the jet power, Q, in RLQs proportional to both the square of the deficit 
in EUV emission, /^, and the accretion rate. The most straightforward interpretation is 
to assume that magnetic flux tubes occupy a fraction, /, of the EUV emitting region since 
this explains the EUV deficit and jet power connection. Furthermore, it was argued that 
the MAA/EMAA variants of this idea explains the precise scaling laws that are implied 
by the observational data. In the following, the experiments designed to test the MAA 
hypothesis are summarized. The speculative nature of the MAA scenario is critiqued as well 
as alternative explanations of the observed scaling laws. 
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In the first experiment, the scatter in Figure 3 indicates that Q/Lboi is strongly corre¬ 
lated with ctEuv as expected from Equation (4), i.e., the jet launching region displaces the 
EUV emitting region in MAA. Secondly, it was shown that the correlation of Q/Fboi with 
oeuv is stronger than the correlation of Q with oeuv- If was argued that this is a statisti¬ 
cally signihcant difference because a partial correlation analysis indicates that the correlation 
of Q/Fboi with ttEuv is statistically signihcant and the correlation of Q with oeuv is spu¬ 
rious. This is evidence of the precise details of MAA that is elucidated by Equation (3), 
ram pressure regulates the strength of the large scale poloidal held strength in the magnetic 
islands. 


In the second experiment. Figure 4 shows a data scatter that is consistent with the jet 
strength being proportional to the square of the EUV dehcit estimated from HST spectra in 
RLQs. The htted slope of the log-log relationships between Q/Lboi and / (Figure 4 is plotted 
for a lower cutoh / = 0.025) are consistent with the prediction of Equation (11) of a slope 
equal to 2. Namely with reasonable variation due to the lower cutoh (0.01 < / < 0.05) the 
htted slope was found to be 2.05 ± 0.17. This is signihcant evidence that the fundamental 
details of the crude MAA model are consistent with the observations. Figure 4 indicates 
that / is a few percent for the RLQs with weaker jets and ~ 40% for the RLQs with the 
most powerful jets. 


It is not claimed that this is the only possible explanation of the EUV dehcit in RLQs. 
Other explanations based on numerical and theoretical models include, lower a in RLQs 


(larger innermost stable orbit), or stronger quenching winds in RLQs per the model of Laor 


and Davis (2014). Note that a scenario of larger black hole mass and lower accretion rates 


in RLQs was ruled out empirically as a plausible explanation due to the indistinguishable 
SED peak in RLQs and RQQs (Punsly 2014). However, none of these other explanations 
naturally produces a radio jet and the correlation in Figure 3. 


Furthermore, it is also not claimed that this is the only possible explanation of Q being 
controlled by Pram- Magnetic hux pinned to the event horizon as in McKinney et ah (2012); 


Tchekhovskoy et al. (2011) or trapped in an equatorial gap between the accretion how and 
the black hole as in Punsly (2008) would also result in this same correlation. However, 
neither of these scenarios explain the Q/Lboi scaling with the square of the EUV dehcit, 
found in Figure 4. A possible reconciliation is provided near the end of this section. 


In Section 2, it was noted that “if the magnetic islands are not extremely short-lived 
transient features in the inner accretion how, an approximate balance of the ram pressure and 
the magnetic pressure of the poloidal magnetic held in the islands must exist.” Currently, the 
dihusion rate of plasma onto and oh of magnetic held lines and magnetic reconnection rates 
are not well known near black holes. Not only are these issues critical for the formation of 
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the magnetic islands, but the time evolution of the magnetic islands is determined primarily 
by diffusion (Igumenshchev 2008 Punsly 2015). These dynamical elements only occur as a 


consequence numerical diffusion in modern simulations in over simplihed ideal MHD single 


fluid models of the physics (Punsly 2015). It is not even clear theoretically what the basic 


principles required for an accurate physical depiction would be. Many issues that are related 


to these topics are active areas of investigation in solar and fusion physics (Bauman et al 2013 


Malakit| 2009 Threlfall et al. 2012 Yamada 2007). As such, it cannot currently be claimed 
or refuted that long-lived magnetic islands can or cannot exist in the inner accretion flow. 
The magnetic islands should become Kruskal- Schwarzschild unstable if the density becomes 
low enough and should move outward, not inward. A key physical element is the numerical 
approximation of the physics that determines the density. If the numerical diffusion is 
signihcant, low enough density may not be achieved before the flux tubes are dragged across 


the event horizon. The simulations in Igumenshchev (2008); Punsly et al. (2009) indicate a 


population of magnetic islands within the innermost accretion flow consistent with the range 
of hlling factors, /, that are found in Figure 4. In the innermost accretion flow, the time 
evolution of the islands is non-steady. The inner accretion flow is a dynamic region. There 
are epochs with very few magnetic islands and epochs in which the (large) magnetic islands 
near the black hole become buoyant after losing mass to the jet and slowly wind their way out 
against the ram pressure of the accretion flow. By contrast, the simulations in [McKinney et 


al. (2012); Tchekhovskoy et al. (2011,2012) that are heavily seeded with large scale magnetic 


flux are devoid of magnetic islands close to the event horizon. This is evidenced by the claim 


m 


McKinney et al. (2012) that no signihcant Poynting hux emerges from this region 


,see 


Equation 2, above) as well as the linked online videos of the simulations. The videos show 
the innermost signihcant, modest, magnetic island concentrations are located at r > lOM 


and they are extremely transient. Thus, even though the simulations in McKinney et al. 


(2012); Tchekhovskoy et al. (2011, 2012) have strong Poynting jets from the event horizon. 


they are not consistent with the EMAA model. This either means that the interpretation of 
the EUV dehcit presented here is wrong or the simulations do not represent the magnetic hux 
evolution accurately. The latter option is favored empirically, since magnetic hux displacing 
EUV emitting gas in the inner accretion how explains two correlated phenomena, jet power 
and the EUV dehcit. It not only explains the correlation, but it predicts the explicit scaling 
laws for jet power with accretion rate and with the degree of EUV suppression implied by 
the observations. The simulations without signihcant vertical magnetic hux in the innermost 
accretion how explain only the scaling law for jet power with accretion rate. 


Recall the fact that the only diherence between the RLQ and RQQ composite continua 
in Figure 1 is the EUV. Comparing this diherence to numerical simulations indicates that 
the magnetic islands are concentrated between the event horizon and an outer boundary 
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of < 2.8M (in geometrized units) for rapidly rotating black holes (jPunsly 2014 Penna et 


al 2010). Thus the simple, homogeneous MAA model does not preclude a black hole spin 


assisted component that displaces the EUV emitting gas inside r < 2M as in the ’’erogspheric 
disk” that has been found to occur in some 3-D numerical simulations of high spin black 
holes (Punsly 2008 Punsly et ah 2009). It is also possible that a certain fraction of the flux 
distribution in the inner accretion flow threads the event horizon and extracts spin energy 
(Blandford and Znajek 1977). In particular, what is shown here is that a significant magnetic 
flux in the inner accretion flow (based on Figure 4 from 2.5% to 45% of the EUV emitting 
region in RLQs is hlled with magnetic islands) explains the EUV dehcit and the source of 
the jet power is proportional to the square of this putative flux. However, the distribution 
of flux can be larger than just the inner accretion flow. For example, the magnetic flux 
distribution might in general thread the horizon and the ergospheric equatorial plane as 
well. This is the EMAA model that was sketched out in the introduction. The magnetic 
islands in the innermost accretion flow would represent the outermost portion of this flux 
distribution in this scenario. This inner region of magnetic flux would be proportional to the 
flux that threads the EUV region for the putative generic flux distribution, oc ^euv- 
The power from the event horizon jet and/or the ergospheric disk jet could be larger than 
the jet power from the EUV region. In this scenario, the flux in the innermost accretion flow 
is merely a “tracer” for flux contained in these interior regions and not coincident with the 
primary source of the Poynting flux that powers the jet. However, the flux permeating the 
innermost accretion flow must be significant and the flux in these interior regions is presumed 
to scale with that in the EUV region. Namely, from Figure 4, for a weak jet like that in 
PKS 1252-1-119, / ~ 4% and for a strong jet like that of 1857-1-566 / ~ 45%. By the EMAA 
assumption, one would expect the event horizon/ergosphere flux to be ~ 10 times larger 
in 1857-1-566 compared to PKS 1252-1-119 on average (and therefore Q/Lboi ~ 100 times 
larger) in order to explain the observational results (Figures 3 and 4) that are described in 
this paper. The implication is still that simulations of jets from these regions would need to 
have a signihcant flux in the innermost accretion flow to be consistent with the observations. 
As with the discussion of the time evolution of magnetic islands above, the existing numerical 
algorithms are single fluid MHD approximations to the physics in which numerical diffusion 
determines reconnection rates and diffusion rates and may not be reliable depictions of the 
physical model of magnetic field dynamics in the innermost accretion flow near black holes. 
As such, a different mathematical description of the reconnection process can lead to very 
different poloidal magnetic flux distributions near the black hole ( Punslyj2015 ). Hence, much 
of this discussion is speculative. On a less speculative note, the fundamental deduction drawn 
from the EUV observations is that future modeling and theory of quasar jet origins should 
contain the feature of an innermost accretion flow threaded by substantial vertical magnetic 
flux. 
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It is necessary for a critical analysis that one segregates the various issues addressed in 
the MAA/EMAA model and the numerical simulations by the degree to which they have 
been verified by the scientific method, observation. The hierarchal list below is in order of 
the degree to which each of these issues is verified by observations. 


1. EUV deficit and radio loudness: The EUV deficit is verified observationally here and in 
Punsly (|2014) to correlate with low frequency radio luminosity on super-galactic scales 


(converted to jet power here in order to make contact with theoretical treatments). 

2. EUV originates from inner accretion flow: Being on the high frequency tail of the opti¬ 
cally thick thermal spectrum this is the obvious interpretation. However, observational 
verification is scientifically much more robust. The only direct method of estimating 
the size of the EUV region is through time variability arguments. This requires the 
ability to collect a signihcant EUV flux on short time scales. The only reasonable set 
of observations related to an active nucleus are the EUVE (Extreme Ultraviolet Ex¬ 
plorer) observations of NGC 5548. Simultaneous UV and EUV monitoring indicated 
that the most likely interpretation of the EUV emission was the Wien tail of the op¬ 


tically thick thermal emission and it had signihcant variability (Marshall et ah 1997). 


Further EUV monitoring found signihcant variability at the smallest sampling time 
scales tyar < 5600 s (Haba et al 2003). Standard arguments based on the light travel 


time across the EUV emitting region and black hole reverberation mass estimates of 
the central black hole indicate the EUV emitting gas is located in a volume with a 


radius, r < 10.5M (Bentz et al. 2007 Denney et al, 2010). The light curve in Haba 


et al (2003) seems to indicate that with a higher sensitivity telescope, the time scale 


for variability would likely be less than 5600 s. This is observational evidence that the 
EUV is radiated from the innermost accretion how. 

3. The only know energy source for the extreme powers in relativistic quasar jets is 
Poynting hux (independent of the details of its origins) and this scales with the square 
of the enclosed poloidal magnetic hux. 

4. Magnetically arrested accretion occurs in some numerical models, however the dy¬ 
namics of the arresting magnetic hux tubes dihers from what is indicated in the 
MAA/EMAA model to varying degrees. The diherence is small to modest in Igu- 


menshchev (|2008); Punsly et al. (2009), that agree with the hll fraction, /, expected 


from Figure 4 in the innermost accretion how. The time evolution generally agrees 
with the MAA/EMAA model except that some of the islands may not agree due to 
some very unstable magnetic outbursts at the smallest radii. The diherence from 


the MAA/EMAA explanation of the observations is large in McKinney et al. (2012); 
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Tchekhovskoy et al. (2011, 2012) for which £11 fraction, / 0 in the innermost accre¬ 


tion flow. The “magnetically choked accretion flows” of McKinney et al. (2012) heat 
the innermost accretion flow compressively as part of the magnetic choking process. 
This likely eliminates the EUV region by drastically elevating the temperature of the 
innermost accretion flow as opposed to suppressing the emission from the innermost 
accretion flow as the observations indicate. Alternatively, depending on the precise, 
unknown, details of radiative transfer in this hot, dense gas, the compressive heating 
might drastically harden the high frequency tail of the optically thick spectrum in 
RLQs, the opposite of what is observed. 


This list is the basis of the logic of the analysis presented in this study. Items 1 and 2 above 
state that observations directly indicate that jets disrupt (suppress) the EUV emission from 
the inner accretion disk and the amount of disruption scales with the power of the jet. The 
third point states that the only known mechanism for driving such a powerful relativistic jet 
is Poynting flux that requires significant poloidal magnetic flux at its source. Thus, with very 
little speculation, it is indicated that both the large scale poloidal magnetic flux at the base 
of the jet and a mechanism that suppresses (but does not eliminate) the EUV emission from 
the innermost accretion flow coexist at the heart of the central engine of RLQs. The main 
assumption of the MAA/EMAA idea is that this is too large of a coincidence, some of this 
magnetic flux must be the same element that disrupts, but does not eliminate the innermost 
accretion flow. If a numerical effort cannot reproduce this circumstance then perhaps the 
numerical approximation to the relevant physical processes requires further development. 
Observations should lead the numerical work. 
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Appendix 1.: Notes on Individual Sources 

0024+224. The NVSS image at 1.4 GHz shows a predominantly core-jet morphology. 
The one-sided jet is likely Doppler boosted kpc emission as is typical of most RLQs viewed 


close to the jet axis Punsly (1995). There are two faint flux density peaks just beyond the hot 


spot at the end of the jet. This is taken as evidence of isotropic halo type emission and this 
flux density is used in the estimate of Q. Even though the radio source is core-dominated, 
the optical polarization and optical variability are small and the emission lines have high 
equivalent widths. Thus, Doppler boosted jet emission in the EUV is considered negligible. 

0232-042. This is a strong lobe dominated source, so the low frequency total flux density 
at 178 MHz is the best estimator for Q. 

0637-752. This is a famous core plus one-sided kpc jet dominated southern hemisphere 
radio source. There is clearly considerable Doppler beaming. The most reliable estimate 
of isotropic flux is to take twice the counter lobe flux density as discussed above. This 
is attained from the ATGA 4.8 GHz radio image Tingay et al (1998). Even though the 


radio source is core-dominated, the optical polarization and variability are small and the 
emission lines have high equivalent widths. Thus, Doppler boosted jet emission in the EUV 
is considered negligible. 

0743-67. This is one of the most intrinsically luminous quasars and is also a very strong 
radio source, both the radio core and the lobes. The lobe flux density is estimated from 


This preprint was prepared with the AAS I+rgN macros v5.2. 
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ATCA radio observations Punsly and Tingay (2005). Even thongh the radio core exceeds 
1 Jy, the optical polarization and variability are small and the emission lines have high 
eqnivalent widths, Thns, Doppler boosted jet emission in the EUV is considered negligible. 
This sonrce has significantly higher visnal extinction than any other sonrce in the sample, 
Ay = 0.362. Thns, nnlike the other sonrces, the detailed form of the Galactic extinction law 
significantly affects the EUV continnnm fit. The average Galactic valne of Ry = 3.1, gives 


a poor fit to a power-law continnnm for the GGM model Gardelli et ah (1989). The reason 
is clear, de-convolving the extinction dne to the A2175 bnmp has created an “artificial” kink 
in the spectrnm presented in this Appendix (this eqnates to 865 A rest frame). A larger 
valne of Ry will remove this artificial bend and these valnes have been fonnd for varions 


lines of sight in the Galaxy Gardelli et ah (1989). Using Ry = 4.5 or Ry = 5.5 makes the 


continnnm look more like a power-law in the spectra above. Even thongh the largest valne 
of Ry gives the smoothest continnnm fit, and improves the scatter plots in Figures 3 and 4, 
a conservative intermediate choice is taken with a large uncertainty assigned in Table 2. 


0959-1-6827. Fortunately, this weak radio source was observed the VLA Landt et ah 


(2006). This is a triple radio source. Such weak sources are very rarely imaged with high 


resolution VLA (as is required at intermediate and high redshift in order to resolve most 
radio sources). 

1022-1-194. This source has a strong core (480 mJy) and a large amount of diffuse, 
elongated emission (160 mJy) at 1.4 GHz in FIRST images. The core spectrum is fiat from 


1.4 GHz to 5 GHz Hutchings et ah (1988). The total spectrum is very steep at low frequency. 


so it seems that the diffuse extended emission dominates at 178 MHz. Thus, the total 178 
MHz flux density is used in the estimate of Q. 

1040-1-123. 3G 245 is the rare quasar with a 1 Jy radio core and very strong symmetric 
radio lobes. In order to get a good estimate of Q, a high resolution radio map is needed 
to extract the radio core and jetted emission Murphy et al. (1993). Even though the radio 


core exceeds 1 Jy, the optical polarization and variability are small and the emission lines 
have large equivalent widths. Thus, Doppler boosted jet emission in the EUV is considered 
negligible. 


1137-1-660. 3G 263 is a lobe dominated quasar, so the low frequency total flux density 
at 151 MHz is the best estimator for Q. 

1229-021. This is a classical triple that is steep spectrum at low frequency. Thus, the 
160 MHz flux density is used to estimate Q. This is a low optical polarization quasar with 
modest optical variability and the emission lines have large equivalent widths. It is concluded 
that the continuum in the EUV is dominated by the optically thick thermal emission. This 
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is corroborated by negligible changes in the EUV continuum between HST observations 
separated by more than 2 years. 

1241+176. A faint, very distant secondary is connected to the radio core by a thin 


bridge in the 5 GHz radio image Kellerman et al. (1994). The flux density of this component 


as determined by the NVSS image at 1.4 GHz is used to estimate Q. HST observations 
separated by 9 years show negligible EUV variability. This combined with large emission 
line equivalent widths indicate that the radio core does not contribute signihcantly to the 
EUV continuum. 

1244+324. This is a classical triple, lobe dominated radio source. Thus, the 151 MHz 
data is used to estimate Q. 

1252+119. This is a core dominated radio source with a weak jet. However, the 1.4 GHz 
radio image shows diffuse emission beyond the end of the jet that is considered evidence of 


a core-halo type source Price et al. (1993). The faint extension is used to estimate the halo 


(lobe) flux density for use in the computation of Q. This is a low optical polarization source 
that is not highly variable in the optical. There is no EUV variability over a 4 year span of 
HST observations. Thus, it is concluded that even though the source is core dominated, the 
relativistic jet contributes negligibly to the EUV. 

1340+289. At 1.4 GHz, the flux is equally split between a “core” and the southern lobe 


Hutchings et al. (1988). The 4.8 GHz image shows a small northern extension, likely part 


of a lobe emission. The overall size is at least 25 kpc at 4.8 GHz making it suitable for this 


sample and the estimation techniques for Q Price et al. (1993). The spectrum is steep at 


low frequency, so it is assumed that the lobe emission dominates the 151 MHz flux density 
and this is used to estimate Q. The optical polarization is low and there is no evidence of 
strong optical variability. Gombined with the large equivalent widths of the emission lines 
implies that there is very little synchrotron contamination of the EUV continuum. 

1340+606. 3G 288.1 is a lobe dominated radio source, so the 151 MHz flux density is 
used to estimate Q. 

1354+19. This source has a strong radio core and a strong one-sided jet, similar to PKS 
0637 -752. Since the counter-lobe is more luminous than the lobe on the jetted side in the 1.4 
GHz high resolution radio images, it is safe to assume that Doppler beaming is insignihcant 


in the jetted lobe Murphy et al.| (1993). Thus, the total lobe flux density at 1.4 GHz (less 
the Doppler beamed one - sided kpc jet emission) is used to compute Q. This object has low 
optical polarization and is mildly variable in the optical. Gomparing the G160L and G270H 
spectra taken 14 months apart, there is no difference in the overlap region, thus the rest 
frame far UV flux is not highly variable. Thus, it is concluded, in spite of the strong radio 
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core, that the EUV is predominantly optically thick thermal emission. 


1415+172. This is a lobe dominated triple Hutchings et ah (1988). Thus, the total low 
frequency flux density is used to estimate Q. 


1857+566. This is a powerful steep spectrum triple Barthel et ah (1988). Thus, the low 
frequency total flux density is used for estimating Q. 


2149+212. This is a compact (20 kpc) steep spectrum triple Barthel et ah (1988). Thus, 
the low frequency total flux density is used for estimating Q. 


2340-036. This source is a symmetric core dominated triple radio source in FIRST 
images. The 1.4 GHz flux density of the lobes is extracted from the FIRST data and is 
used in the computation of Q. Fven though the radio source is core-dominated., the optical 
polarization and variability are small and the emission lines have high equivalent width. 
Thus, Doppler boosted jet emission in the FUV is considered negligible. 


Appendix 2.: HST Spectra 

Figure 5 are the FUV spectra corrected for Galactic extinction, the Lyman a forest and 
Lyman limit systems as discussed in the Sample Selection section. The spectra are arranged 
in order of increasing Q/Lboi in order to see the trend of increasing FUV deficit regardless 
of the precise continuum power law fits (the black lines). The spectra are log-log plots of 
Fx as a function of A. The data was downloaded from MAST and smoothed as required to 
enhance the dehnition of the continuum 
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Fig. 5.— a 
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Fig. 5.— b 
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Fig. 5.— c 
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Fig. 5.— d 
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Fig. 5.— e 
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Fig. 5.— f 
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3C 263 
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1022+194 



Fig. 5.— g. The spectral data for 3C 263 was generously provided by Michael Shull and 
Matt Stevans 
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Fig. 5.— h. 
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Fig. 5.— i. 
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Fig. 5 
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Fig. 5.— k. 





























































































































